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. In the framework of the Schwinger-Dyson equation and the Bethe-Salpeter 

equation in the improved ladder approximation, we calculate the S parameter 
and an analogue of the tt + — tt° mass difference AmJ = m^. + — m^ as well as 
^ ■ the NG boson decay constant f n on the same footing in the large Nf QCD, 

through the difference between the vector current correlator Hyv an d the axial- 
vector current correlator Haa- Approaching the chiral phase transition point 
a* — ► a C r(= 7r /4) from the broken phase, where a* is the the gauge coupling 
on the infrared fixed point, AraJ as well as goes to zero with the essential- 
singularity scaling (Miransky scaling), while the ratio indicates a blowing up 
enhancement reflecting the characteristic behavior of the large Nf QCD as a 
walking theory which is expected to scale as Am 2 // 2 ~ (a*/a cr — l) -1 ^ 2 . On 
the other hand, the S parameter takes values somewhat smaller than that of 
the real-life QCD and indicates slightly decreasing tendency as we approach 
the phase transition point. 



*A preliminary report was given at the 2004 International Workshop on Dynamical Symmetry 
Breaking (DSB04), December 21-22, 2004, Nagoya Univeristy, Nagoya. 
t Present address 



1 



1 Introduction 



The electroweak symmetry breaking in the standard model (SM) is well described by 
introducing the Higgs boson. However, the Higgs sector in the SM cannot explain the 
origin of the electroweak symmetry breaking, which means the SM cannot explain 
the origin of mass. Moreover, the mass of the Higgs particle has to be fine tuned 
since it receives quantum correction which is proportional to the square of the cutoff 
scale A. These facts naturally lead us to the idea that the Higgs sector in the SM is 
nothing but an low energy effective theory of some underlying theory. To consider 
the existence of a new strong dynamics is one example of this idea. 

The "large Nf QCD", a jargon of the 577(3) gauge theory with a large number 
of massless fermions (N f < fN c ) El S 13 13 13 13 HOJ, is one of the candidates 
for such a strong dynamics like the technicolor ^T] which breaks the electroweak 
symmetry. Actually, it has interesting features: 

• A dynamical model for the walking/conformal (scale-invariant) technicolor 

It was observed [2] that QCD has an infrared (IR) fixed point a = a* for a 
large number of massless fermions (Nf < y^ c ) a ^ ^ ne t w °T°°P beta function. 
It was regarded [3 El as a good example of the walking technicolor jT^l , because 
the running gauge coupling near the infrared fixed point is almost constant, 
a(Q 2 ) ~ a*, namely "walking" for a wide range of the momenta Q. 

Moreover, it was found [3| through the improved ladder Schwinger- Dyson (SD) 
equation that chiral symmetry restoration takes place for certain Nf such that 
the IR fixed point becomes less than the critical coupling which is determined 
by the SD equation, namely A^ rit < N f < fiV c , where A^ rit ~ AN C (= 12 for 
N c = 3). In Ref. JT] this chiral phase transition at Nf 11 was further identified 
with the "conformal phase transition" which was characterized by the essen- 
tial singularity scaling (Miransky scaling). Moreover, the chiral restoration in 
large Nf QCD was also observed in lattice simulations [H] and by other various 
methods such as the dispersion relation instanton calculus (Zj, effective field 
theoretical approach jS], renormalization group flow equations [TU], perturba- 
tive calculus , etc., although the value of A^ nt is raging from Nf Tlt ~ 5 to 
Nf nt ~ 12 (N c = 3) depending on the approaches. 

• A model for the electroweak baryogenesis 

The electroweak baryogenesis [15] (for a review see, e.g., Ref. [Tfi] ) is an attractive 
scenario but has serious difficulties: 1) lack of mechanism to produce the first 
order phase transition consistently with the present experimental lower bound 
of the Higgs boson mass and 2) smallness of the baryon asymmetry if the CP 
violation is solely due to the KM phase in the Standard Model. The large 
Nf QCD when applied to the technicolor has new ingredients to solve these 
problems: It is widely believed [TJj that QCD phase transition is the first order 
when the number of massless flavors exceeds 3 and, we expect, so is the chiral 
phase transition in the large Nf QCD. Moreover, it was argued [IB] that the 
technicolor may have extra sources of the CP violation other than the KM phase 
and supply enough amount of CP violation needed for the baryogenesis. 
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Besides application to the technicolor, much attention has been paid to the prop- 
erty of the phase transition not just the existence of the chiral phase transition in 
the large Nf QCD. Especially, it is interesting to ask what are the light degrees of 
freedom near the phase transition point in the large Nf QCD: For example, in the 
manifestation of the chiral symmetry restoration a la linear sigma model (Ginzburg- 
Landau type effective theory), the scalar bound state becomes a chiral partner of the 
pseudoscalar bound state and becomes massless at the phase transition point when 
approached both from the broken and the symmetric phases. However, it was em- 
phasized from the viewpoint of the conformal phase transition ||] that the Ginzburg- 
Landau type effective theory breaks down in the approach of the improved ladder SD 
equation, signaled by the absence of massless scalar bound state in the symmetric 
phase. Quite recently, such a peculiarity of the phase transition was also observed 
based on the renormalization-group analysis ^U]. Besides, from the viewpoint of the 
conformal phase transition , it is natural to suppose that all the existing bound states 
become massless near the phase transition point when approached from the broken 
phase (see Ref. |19j). On the other hand, in the vector manifestation [§| obtained 
by the effective field theoretical approach based on the hidden local symmetry the- 
ory pOj, it is the vector bound state which becomes massless as a chiral partner of 
the pseudoscalar bound state. It remains unclear at this moment whether or not the 
vector manifestation contradicts the conformal phase transition. 

In Ref. [21], in the framework of the SD and the homogeneous Bethe-Salpeter 
(HBS) equations (see [22] for a review), we studied which types of the light bound 
states actually exist near the phase transition point, and investigated the critical 
behavior of their masses directly from QCD. We found that approaching the chiral 
phase transition point from the broken phase, the scalar, vector, and axial-vector 
meson masses vanish to zero with the same scaling behavior, all degenerate with the 
massless pseudoscalar meson. 

In this paper, we further investigate the properties of the chiral phase transition 
in the large Nf QCD through the critical behavior of the S parameter and the "7r + — 
7T° mass difference AmJ = — m^ ", by extending our previous work [22] on 
the S parameter and AmJ done for the real-life QCD with N c = Nf = 3. They 
are interesting quantities when we apply the dynamical symmetry breaking scenario 
based on the large Nf QCD for the model buildings beyond the SM like the walking 
technicolor, UV completions (underlying theories) for the little Higgs models [25 an d 
the Higgsless models j2HJ, etc. 

Actually, there are strong constraints on the electroweak S parameter [26j from 
electroweak precision tests, it is quite important to estimate the S parameter theoret- 
ically when we consider the technicolor-like scenario. Estimation of masses of pseudo 
NG bosons is also important to check that the model is consistent with the fact that 
such a pseudo NG bosons have not been discovered so far. Then, it is desirable to 
estimate the S parameter and Am^ reliably. 

Several attempts to estimate the parameter S, which is the contribution to the S 
parameter from one weak fermionic doublet, for the walking technicolor theories and 
the large Nf QCD have been made so far by using several methods [2Z1I2H]- However, 
it is difficult to estimate the accuracy of these calculation (see, for example, Ref. 29J). 



3 



What makes it difficult to estimate the S parameter is its strong dependence on the 
nonperturbative dynamics of the low momentum region. Then it is quite important to 
calculate the S parameter by using the nonperturbative method in a way to directly 
deal the fermions and the gauge bosons as the fundamental degrees of freedom. 

As such, in the previous work , we calculated the S parameter and AmJ in the 
real-life (Nf = 3) QCD through the difference between the vector current correlator 
Hvv an d the axial- vector current correlator 11^. Hyv ~ ^-aa were calculated in the 
framework of the SD equation and the inhomogeneous Bethe-Salpeter (IBS) equation 
in the improved ladder approximation. It was stressed that these physical quantities 
are calculated by using the correlators in the space-like region, so that no analytic 
continuation is needed. The results in Ref. shows that both the S parameter 
and Am^ can be fit to the experimental values at the same time. By fitting to the 
calculated data using the pole saturated form of Hyv ~ ^-aa, we also derived masses 
and decay constants of p meson and A\ meson, which we found are consistent with 
the experiments. This is quite encouraging to extend our method to the large Nf 
QCD. 

In extending the previous method to the large Nf QCD, let us consider the sit- 
uation that SU(Nf) L ® SU(Nf) R chiral symmetry is explicitly broken by the U(l) 
gauge symmetry. We take the generator of this U(l) gauge symmetry as 



Q 



/2/3 



V 



-1/3 



/ 



Then, we calculate the mass of the pseudo Nambu-Goldstone (NG) boson which is 
associated with the following generator: |5U| 



T, 



pNGB 
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/ 
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In the case of the real-life QCD, the above pseudo NG boson corresponds to the n + 
meson, and its (squared) mass caused by the U(l) electro-magnetic gauge interaction 
corresponds to the tt + — 7r° mass difference AmJ = — m^ . Then, throughout 
this paper, we call the squared mass of the above pseudo NG boson "AmJ" even in 
the case of Nf ^ 3. 

We then calculate S and AmJ as well as f n at the same time. We find that as 
we approach the critical point both Am^ and f% vanish in the essential singularity 
scaling (Miransky-type scaling), ~ exp (—2q/ \J a^/a cr — , (q = const.) with the 
ratio AmJ//j being much enhanced compared with that of the real-life QCD and 
having a blowing up tendency reflecting the inverse square root scaling, AmJ//^ ~ 



V 



-1/2 



which is characteristic to the walking gauge coupling with large 



anomalous dimension 7 m ~ 1 of the large Nf QCD. In fact the integral for the AmJ//^ 
is expected to be logarithmically divergent for the walking theory, ~ In (A 2 /A^J, if 
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the walking behavior Hyv — ^-aa ~ 1/Q 2 is operative for A x < Q 2 < A 2 , where A is 
the (two-loop) scale parameter of the large Nf QCD like Aqcd i n the real-life QCD, 



and A x ~ 4:7if n ~ Aexp ^—q/ya^/a,,-,. — 1J. In the actual numerical calculations we 

used an artificial cutoff A^ um which is defined for the sake of numerical calculations 
as Hvv — II,4a|q2 = a 2 = ^ \^vv — Ha a] |q 2 =o ( we will later include contributions 
from A^ um < Q 2 < A 2 , based on the analytical considerations on the walking theory). 
Up to this limitation the numerical result is already enhanced up till Am 2 // 2 ~ 0.6, 
about four times larger than that obtained in the same method for the real-life QCD, 
Am 2 // 2 ~ 0.123 |23! (experimental value: Am 2 // 2 = 0.148 ± 0.001), even though 
our calculation is still away from the critical point for the reason of the computation 
ability. (If we include contributions from A^ um < Q 2 < A 2 , the enhancement is more 
dramatic) . 

This implies a large enhancement factor for the pseudo NG boson mass in the 
walking/conformal (scale-invariant) technicolor [13] • 1 If it is applied to the popular 
one-family model (Farhi-Susskind model) |12j . the mass 2 of the P ± (and also of the 
colored pseudo NG bosons) as a scale-up of Am 2 essentially irrelevant to the ETC 
interactions, will be enhanced also by a large anomalous dimension 7 m ~ 1 by a factor 
four (to nine): The mass of will be boosted to typically more than 200 (— 300) GeV 
mass range instead of lOOGeV the value in a simple scale up of QCD. If we consider 
the model which is closer to the critical point, the enhancement factor should be more 
dramatic due to the inverse square root scaling mentioned above. 

On the other hand, the values of S in our calculations are somewhat smaller 
S = 0.25 — 0.30 than that obtained in the same method for the real-life QCD (N c = 
Nf = 3), S = S ~ 0.33 [23] (experimental value: S = 0.32 ± 0.04), and indicate a 
tendency to decrease slightly as we approach to the critical point. Although it is still 
larger than the simple one-loop result S = ^ = 0.16 (N c = 3), it could be suggestive 
for a further decreasing towards the critical point. It would motivate further study in 
the points closer to the critical point to see whether or not technicolor model based on 
the large Nf QCD makes S consistent with the present experimental bound S < 0.1. 

The calculation in this paper turns out to be the first example which estimates 
the S parameter and Am 2 in the large Nf QCD as a model of walking technicolor by 
directly solving its dynamics 2 . 

This paper is organized as follows. In the next section, we briefly review the large 
Nf QCD. In section [3] we introduce sum rules: The DMO sum rule or the zeroth 
Weinberg sum rule for the parameter S, the first Weinberg sum rule for / 2 , and the 
DGMLY sum rule or the third Weinberg sum rule for Am 2 . We then rewrite them in 
terms of the current correlators Hvv — ^-aa- In section 0] we show how the current 
correlators are obtained from the BS amplitudes calculated from the IBS equation 
in the space-like region. In section we introduce the IBS equation. In section El 
we show the results of numerical calculation for critical behaviors of /„-, Am 2 and S. 

X A similar enhancement of pseudo NG boson mass in the walking technicolor (other than the 
large Nf QCD) was also observed by a numerical study in a different approach 

2 In Ref. [32], the S parameter in the walking technicolor theory was calculated by using the BS 
equation. However, the S parameter in the large Nf QCD was not calculated. 
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Section [7| is devoted to conclusion and discussions. 



2 Large N f QCD 

In QCD with Nf flavors of massless quarks, the renormalization group equation 

(RGE) for the running coupling a(/x) ( = ^-j^- J in the two-loop approximation is 
given by 

/i— = —ba 2 (/j,) — ca 3 (/i), (2.1) 
dfx 

where 

1 1 / N 2 -1 \ 

b = - (11N C — 2Nf) , c = — - 10N c N f - 3^— N f j . (2.2) 

From the above beta function we can easily see that, when b > and c < 0, i.e., Nf 
takes a value in the range of Nf < N f < ^N c {N* f ~ 8.05 for N c = 3), the theory 
is asymptotically free and the beta function has a zero, corresponding to an infrared 
stable fixed point 0111], at 

«* = - — . (2.3) 

c 

Existence of the infrared fixed point implies that the running coupling takes a 
finite value even in the low energy region. Actually, the solution of the two loop RGE 
in Eq. (J2.1j) can be explicitly written [331 EH! * n a ^ the energy region as 

a (jj) = a, [ W{n ba */eA ba *) + 1 ] _1 , (2.4) 

where W(x) = F~ 1 (x) with F(x) = xe x is the Lambert W function, and A is a 
renormalization group invariant scale defined by jHj 



A = /i exp 



1 , fa* — a([i) 
log 



6 a* y a(/j) J b a(n) 



(2.5) 



We note that, in the present analysis, we fix A to compare the theories with a 
different number of flavors, and that we have no adjustable parameters in the running 
coupling in Eq. (J2.4)) . accordingly. 3 We show an example of a(/i) for Nf = 9 in Fig.^ 
The fact that the running coupling is expressed by a certain function as in Eq. ()2.4j) 
implies that, in the case of the large Nf QCD, we do not need to introduce any infrared 
regularizations such as those adopted in Ref. [23], or in Refs. jSS] for studying real-life 
QCD with small Nf in which the infrared regularization parameter must be chosen in 

3 Note that we do not have to specify the value of A, since in this paper we calculate only 
dimensionless quantities, S parameter and Am^//^ which are independent of the actual value of 
A in unit of, say GeV. When we apply the dynamics to the model building for the electroweak 
symmetry breaking, we should specify the value of f n (of order of the weak scale 250 GeV but 
actually model-dependent) and then determine the scale of A as a function of N c (= 3) and Nf as 
well as (See FigOJ. 
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Figure 1: Two-loop running coupling of the large Nf QCD for Nf = 9. 

such a way that the running coupling in the infrared region becomes larger than the 
critical value a cr = it/A for realizing the dynamical chiral symmetry breaking. The 
running coupling in the large Nf QCD takes a certain value in the infrared region 
for given Nf, so that we can definitely determine, within the framework of the SD 
equation, whether or not the dynamical chiral symmetry breaking of SU(Nf)i, x 
SU(Nf) R — > SU(Nf)y is realized. Actually, the value of a* decreases monotonically 
with increasing Nf, and the chiral symmetry restores when Nf becomes large enough. 
In Refs. jH], it was shown that the phase transition occurs at Nf Tlt ~ 11.9 for N c = 3 
(corresponding to a* = a cr = 7r/4). However, the value of Nf Tlt itself should not be 
taken seriously, since a CT = 7r/4 to be equated to a* is already so large as to invalidate 
the perturbation which determined the value of a*. There is also some ambiguity of 
the ladder approximation which determined the value a cr = tt/A itself. In fact lattice 
simulations [5 j suggest smaller value 6 < Nf nt < 7 and other approaches also suggest 
different values CD] • 

3 Expressions of f n , S and Am^ in terms of current 
correlators 

In this section we introduce sum rules for expressing the order parameter f n , the 
parameter S which is the contribution to the S parameter from one weak fermionic 
doublet, and the square mass of the pseudo NG Am^ in terms of Hvv — ^-aa in the 
space-like momentum region. 

Let us begin with introducing the vector and the axial-vector currents as 

V£(x) = ^(x)T a ^(x) , 

A;(x) = ij(x)T a W ij(x) , (3.1) 

where T a (a = 1, 2, . . . , Nj — 1) is the generator of SU(Nf) normalized as tr(T a T b ) = 
\5 ab . As we discussed in the previous section, for Nf < Nf 11 there occurs the spon- 
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taneous symmetry breaking of SU(Nf\ x SU(Nf)-& — > SU(Nf)y. As a result, the 
massless NG boson tc appears, whose decay constant f n is defined as 



(0\A a (0)\7r b (q))= iq JJ ab , (3.2) 



where a, b = 1, 2, ... , Nj — 1. As is well-known as the first Weinberg sum rule 
the above decay constant is related the difference of the vector current correlator Hyv 
and the axial- vector current correlator 11^ as 

fl = T1 VV (0) - Tl AA (0) , (3.3) 

where the correlators are defined from the currents in Eq. (J3.1|) as 

f* (ML - g \ njJ (g 2 ) = i [ d'x (0\Tj;(x)J h v (0)\0) , (3.4) 



(•£(*) = V£(*MS(*)) . 



When the large Nf QCD is considered as the underlying theory for describing 
the Higgs sector of the standard model as in the technicolor-like scenario , it 

is important to study the S parameter which gets a strong constraint from the 
electroweak precision tests j^Zj. Since entire contributions to the S parameter depend 
on the detail of the model structure, we concentrate on calculating the contribution 
from one S77(2)l doublet with the hypercharge Y = 1/6 and two singlets with Y = 
2/3 and Y = —1/3. This parameter S is related to the above correlators as 



S= -4vr 



d 



dQ 2 



Tl vv (Q 2 )-Tl AA (Q' 



(3.5) 



which is nothing but the DMO sum rule [SB] or often called the "zeroth Weinberg 
sum rule" . 

Mass of the pseudo NG boson is also an important quantity in the technicolor-like 
scenario. In the present analysis, we consider the situation that SU(Nf) L <& SU(Nf) R 
chiral symmetry is explicitly broken by the U(l) gauge symmetry, the generator of 
which is given as in Eq. (|1.1|) . and calculate the the mass of the pseudo NG boson 
associated with the generator given in Eq. (jl.2|) . The squared mass of the above 
pseudo NG boson Am 2 is related to the above correlators as 



Ami = P%r d ® 2 \ U w(Q 2 ) ~ KAAiQ 2 )} ■ (3.6) 
Airf 2 Jo L 

where a em (= |^ = y|y). Equation ()3.6|) was first derived by Das et al. [10] m the 
form written in terms of the spectral function instead of the current correlators and 
may be called the DGMLY sum rule or the "third Weinberg sum rule" . 

As was studied in, e.g., Ref. AmJ is interesting quantity to study because 
it is related to the structure of the vacuum. When we switch off the above gauge 
interaction, we have Nj — 1 massless NG bosons associated with the spontaneous 
breaking of SU (Nf)^ ® SU (A^-)r chiral symmetry down to SU (Nf)y symmetry. The 
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existence of the gauge interaction explicitly breaks the SU(Nf)-^ Cg> SU(Nf)n chiral 
symmetry, which makes some of NG bosons (say tt + ) be massive while others (say 
7T°) remain massless. The interesting point here is whether AmJ becomes positive 
or negative, which is called "vacuum alignment problem" j3T]. Negative AmJ means 
that fluctuation of 7r + field around (tt + ) = is unstable and the vacuum with (7r + ) = 
is not a true vacuum. If this is the case, 7r + has non-zero vacuum expectation value 
and U(l) gauge symmetry is broken. In the case of the real-life (Nf = 3) QCD, we 
have shown 23 i that positive AmJ is successfully reproduced in the framework of the 
SD equation and the IBS equation with the improved ladder approximation. In our 
best knowledge, however, no one has explicitly shown whether positive Am^ is also 
realized in the case of the large Nf QCD, especially near the critical point, so our 
calculation is the first one which explicitly investigates the stability of the vacuum in 
the large Nf QCD. 

4 Current correlators from BS amplitudes 

In the previous section we have written down the sum rules for /„., S and AmJ in 
terms of the current correlators. In this section, we show how the current correlators 
are obtained from BS amplitudes which will be calculated from the IBS equation. 

For deriving several properties of hadrons as boundstates, we often need to perform 
calculations in the time-like momentum region. However, it is difficult to solve the BS 
equation and the SD equation in the time-like region since we have to carry out the 
analytic continuation of the running coupling from the space-like region to the time- 
like region. 4 In the present analysis, on the other hand, we need the BS amplitude 
only in the space-like region in order to calculate the current correlators appearing in 
the sum rules for Am^, f n and the QCD S parameter. 

Now, the BS amplitudes (J = V ■, A) are defined in terms of the three-point 
vertex function as follows: 

«UT")} j -^e- ipr X$(p-,q,e) = e>*J #x<*°(0\T ^/(r/2) (-r/2) Jftz) |0>, 

(4.1) 

where g M is the total momentum of the fermion and the anti-fermion, and p M is 
the relative one. & is the polarization vector defined bye-g = 0, e • e = — 1, 
and (/, /'), (cx.,0) are flavor, color and spinor indices, respectively. Closing the 
fermion legs of the above three-point vertex function and taking the limit r — > 0, we 

4 In general, one-loop running coupling of QCD is not analytic, so that several models for running 
coupling are introduced to solve the SD and BS equations in the time-like region. On the other hand, 
the running coupling at two-loop level becomes fully analytic near the phase transition point in the 
large Nf QCD. In Ref. |21| . we approximated the analytic running coupling by a simple function 

- = a* #(A 2 — (x + y)), and studied the critical behaviors of f n as well as masses and decay 

constants of scalar, vector, and axial-vector mesons. We note that, in the present analysis, we do 
not use the approximate form, but use the fully analytic two-loop running coupling. 
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can express the current correlator in terms of the BS amplitude as follows: 



*W) = I E / ^ f tr [(e • G(-O) X (-0(p; g, e)] (4.2) 

where 

= 7„ G<*> = 7m75 , (4.3) 

and N c = 3 is the number of colors. In the above expression we averaged over the 
polarizations so that Ujj(q 2 ) does not depend on the polarization. 

We expand the BS amplitude Xafl ip'i 9) e ) i n terms of the bispinor bases and 
the invariant amplitudes Xi^ as 

x (J) feg,e)l = Efr^fa&e)] ,xS J W, (4.4) 

J otp — ■ L J ap 

1=1 

where g M = q^j 'y/Q 2 with Q 2 = —q 2 . The bispinor bases can be chosen in a way 
that they have the same properties of spin, parity and charge conjugation as the 
corresponding current J£(x) has. We adopt the following bispinor bases for the vector 
vertex: 

rf° = ^ r 2 ^ = ^](p.$), r^ = ^,|], rf> = ±[^,|] (4.5) 

rr } = (e-p), rj^tfcp), T^ = i(p.q)(e.p), r 8 ^ = ^,|](c P ), 
where [a, 6, c] = a[&, c] + 6[c, a] + c[a, 6]. For the axial- vector vertex we use 

if> = ^75, r^ = ^ 7 5, if> = i[6f](p.$) 7 6, 

ri A) = ^,fl75, r^ = (c.p)(p.g) 7 6, r^=^( e -p) 75 , 

T 7 A) = |(e-p) (p- g) 7 5, ri- A) = ^,|](e.p)(p.g)76. (4-6) 

From the above choice of the bases, we can easily show that all the invariant ampli- 
tudes x\ J ^ are the even functions of (p ■ q) using the charge conjugation property of 
the current. 

In the present analysis we fix the frame of reference in such a way that only the 
zero component of the total momentum q^ becomes non-zero. Furthermore, we study 
the case where g M is in the space-like region. Then, it is convenient to parameterize 
the total momentum q^ as 

g" = (iQ, 0,0,0). (4.7) 

For the relative momentum p M , we perform the Wick rotation, and parameterize it 
by the real variables u and x as 

p ■ q — —Q u , p 2 = —u 2 — x 2 . (4.8) 
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Consequently, the invariant amplitudes xi^ become functions in u and x: 



Xi J) = Xi J \%x;Q). 



(4.9) 



From the charge conjugation properties for the BS amplitude x^ an d the bispinor 
bases defined above, the invariant amplitudes xf\ u ^ x ) are shown to satisfy 



Xi J \u,x;Q) = Xi J \-u,x;Q) . 
Using this property of the invariant amplitudes, we rewrite Eq. (|4.2j) as 



poo poo 

H-vv(Q 2 ) = —t du dx 

7T d JO JO 



X 



J\f poo poo 

Raa(Q 2 ) = —t du dx x 2 

7T JO JO 



-X { P (u, x; Q) + ^-Xe K x ; Q) 



Xi A \u,x;Q) - %rXft\% x iQ) 



(4.10) 



(4.11) 



(4.12) 



Here, we used the expanded form of the BS amplitude shown in Eq. (j4.4j) and carried 
out the three dimensional angle integration. 

From Eqs. (|4.11j) and ()4.12|) . the quantity Hyv — ^-aa is expressed as 



Il vv -U 



AA 



E 



d A p N c 



i{2n) 



tr 



— r / du 



7T 3 JO 



„2 



ix [J \p\ q^)-Hx {A \p\ q,e) 

X ( p (u, x; Q) + xS A) (u, x; Q)) 



+ - 



x 



X* (u, x; Q) + xi A) {u, x; Q)) 



. (4.13) 



We note that, although either Hyv or ^-aa is logarithmically divergent quantity, the 
difference Hyv ~ ^-aa becomes finite due to the cancellation of the divergence ensured 
by the chiral symmetry. 



5 Inhomogeneous Bethe-Salpeter equation 

In this section we introduce the inhomogeneous Bethe-Salpeter (IBS) equation from 
which we calculate the BS amplitude defined in the previous section. 

The IBS equation is the self-consistent equation for the BS amplitude x^ J \ an d it 
is expressed as (see Fig. E] for graphical expression) 

T(p-q)x {J \p;q,e) = e ■ + K(p;k) * X (J \k;q,e). (5.1) 

The kinetic part T is given by 

T{p; q) = Sr\p + q/2) ® S F \p - q/2) , (5.2) 
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Figure 2: A graphical expression of the IBS equation in the (improved) ladder approx- 
imation. 



where Sf is the full fermion propagator iSp 1 (p) = p 1 — S(p). (Note that wave function 
renormalization factor A(p) becomes unity when we adopt the Landau gauge.) The 
BS kernel K in the improved ladder approximation is expressed as 

where g(p, k) is the running coupling of QCD whose explicit form will be shown later. 
In the above expressions we used the tensor product notation 

{A®B) X = A X B , (5.4) 

and the inner product notation 

K(p; k) * X {J) (k; q, e) = j K(p, k) X (k; q) . (5.5) 

The mass function of the quark propagator is obtained from the SD equation (see 
appendix IXj) : 

E{p) = K(p,k)*iS F (p). (5.6) 

It should be stressed that we must use the same kernel K(p, k) as that used in the IBS 
equation for consistency with the chiral symmetry ^21 HEj- Numerical method for 
solving the SD equation and the IBS equation are shown in appendix |B] and appendix 
O respectively. 



6 Numerical results 

In this section, we show the results of the calculations for /„., S, and AmJ in the large 
Nf QCD near the chiral phase transition point. The procedure of the calculation is 
as follows. First, we solve the IBS equation and the SD equation simultaneously, 
and obtain BS amplitude x V ' A - (Details of numerical methods to solve the SD and 
the IBS equations are shown in appendix [B] and appendix 0) Then, by closing the 
fermion legs of BS amplitudes, we obtain current correlators YlyviQ 2 ) an d Haa(Q 2 )- 
(See sectionHJ) Once we obtain U VV (Q 2 ) and Uaa(Q 2 ), we calculate S, and Ami 
from the relations in Eqs. ()3.3|) . ()3.5|) and ()3.6|) . 



12 



fn 
A 



0.001 
0.0008 
0.0006 






























































0090° 







OCcr 



J 0.8 

(11.88) 



0.85 
(11.75) 



0.9 
(11.63) 



0.95 
(11.52) 



1 a* 

(11.42) (Nf) 



Figure 3: Values of /„. calculated from Eq. (|3.3() for several values of a* in the range 
of a* € [0.89 : 1] (indicated by O). Values of Nf corresponding to each a* are shown 
in parenthesis. 



6.1 Critical behavior of f n in the large Nf QCD 

In Fig. El we plot f n for several values of a* in the range of a* G [0.89 : 1]. Note that 
we use a* instead of Nf as an input parameter, because once we choose a value of Nf, 
the value of a* is uniquely determined from Eq. (|2.3J1 . a* = a cr implies Nf = 11.91. 
In Fig. El we have shown corresponding values of Nf in the parenthesis. From this 
figure, we can see that f n goes to zero when we reduce the value of a* to its critical 
value a cr ~ it/ 4 (or increase the value of Nf to its critical value Nf 11 ~ 4iV c ), which 
is consistent with the chiral phase transition described in sectional Furthermore, the 
scaling behavior is identified as essential singularity type, 

U=pAexp\- rJ—^\ (6-1) 

with p = 15.1 and q = 4.9 as the best fit. This is similar to the case with constant 
coupling with sharp cutoff at A [2*T]: 5 

U = dAexp\- ; J r _ i 1 , (d~1.5). (6.2) 



5 At first glance there is some numerical difference between the case with a constant coupling and 
that with a two-loop coupling, which is mainly due to the different meanings of A for these two 
cases. Actually, such a superficial difference also takes place for m = S(m 2 ) which is essentially 
proportional to /„■, with the ratios being the same in both cases. 



13 



2 



A 



2 



7e-07 - 
6e-07 - 
5e-07 - 
4e-07 - 
3e-07 - 
2e-07 - 
le-07 - 




oc cr 



J 0.8 
(11.88) 



0.85 
(11.75) 



0.9 
(11.63) 



0.95 
(11.52) 



i a* 

(11.42) (Nf) 



Figure 4: Values of AmJ calculated from Eq. I|3.6|l for several values of a* in the range 
of a* S [0.89 : 1] (indicated by O). 

6.2 Critical behavior of Am^ in the large Nf QCD 

In Fig.Hl we plot AmJ calculated for several values of a* in the range of a* G [0.89 : 1]. 
In our numerical calculation of AmJ, we introduced the numerical cutoff A^ um : 

Hv-a(Q 2 = A|um) J_ /„ o\ 

Hv-a(Q 2 = 0) 50 1 ' ] 



for the integration of Uv-a(Q 2 

Ami -- 



[n vv - TIaa] (Q 2 ) in Eq. as follows: 



3a, 



r A 2 

I num o , r 

/ dQ 2 U v _ A (Q^ 
Jo 



4tt/2 J q 



(6.4) 



We shall discuss later inclusion of the contributions from A^ um < Q 2 < A 2 . 

From this figure, we see that Am 2 vanishes to zero when approaching the chiral 
phase transition point from the broken phase. To see the scaling behavior of Am 2 
near the phase transition point, we plot the values of Am 2 // 2 in Fig. El for several 
values of a* in the range of a* G [0.89 : 1]. From this figure, we can see that Am 2 // 2 
in the large Nf QCD near the phase transition point gradually increases 

Att? 2 

—^ = 0.4 to 0.6, (6.5) 

J IT 

in the range of a* G [0.89 : 1] investigated here. 

The most outstanding feature of the value in Eq. (J6.5j) is that it is rather large 
compared with that of the real-life QCD 

Ami 

0.123, (6.6) 



f 2 



real-lifeQCD 
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OC C r _f 0.8 0.85 0.9 0.95 1 0C* 

(11.88) (11.75) (11.63) (11.52) (11.42) (Nf) 

Figure 5: Values of AmJ//, for several values of a* in the range of a* G [0.89 : 1] 
(indicated by O). Solid line shows the fitting function /(a*) given in Ea. i|tj.ll[l . Value 
of Am%/f% of the real-life QCD calculated in Ref. is also shown in the figure. 
Here, we used the data with the infrared cutoff parameter of the running coupling 
t F = 0.04. 

calculated in the same method with the numerical cutoff AN um defined as Eq. ()6.3|) 
(The experimental value: Am 2 /f 2 = 0.148 ± 0.001). We shall argue that this en- 
hancement reflects the large anomalous dimension 7 m ~ 1 characteristic to the walk- 
ing/conformal (scale-invariant) technicolor [T3] . 

6.3 Inverse square root scaling of Am^ 

To see the cause of this enhancement of Am 2 /f 2 let us compare the behavior of 
U V _ A (Q 2 )/ fl versus Q 2 /f 2 of the large N f QCD with that of the real-life QCD, 
see Fig. El It is clear from Fig. El that Uv-a{Q 2 )/ f% i n the large Nf QCD is more 
slowly damping in high Q 2 than that in the real-life QCD, which consequently yields 
bigger area integral in Fig. El than that of the real-life QCD. This is the cause of the 
enhancement of Am 2 / f 2 . 

To understand the slow damping it is convenient to divide the integral region 
of the DGMLY sum rule ("third Weinberg sum rule"), Eq. (|3.6|) . into high and low 
energy regions as 

Am * - ii(r <i<?2+ /? <32 ) n "-« 2) 

= Am l(IR) + Am l(VV), ( 6 - 7 ) 

where A x was introduced in such a way that, in the high energy region Q 2 > A^, 
the behavior of [Ilyy(Q 2 ) — ll^v^Q 2 )] is well approximated by the one obtained from 
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Figure 6: Comparison of H\/=a(Q 2 ) between the large Nf QCD and the real-life QCD. 
Solid line indicates Hv-a(Q 2 ) in the large Nt QCD for a* = 0.90, while dotted line 
indicates Hv-a(Q 2 ) in the real-life QCD (Nf = 3). In this figure, both horizontal and 
vertical axes are normalized by the respective value of f 2 . Hv-a(Q 2 ) in the real-life 
QCD was calculated in Ref. |2Sj- Here, we used the data with the infrared cutoff 
parameter of the running coupling tp — 0.04, which reproduce the QCD S parameter 
and Am^ in good agreement with experiments |23j. 
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the operator product expansion (OPE) technique in the high energy region (see, e.g., 
Ref. jSl). In the real-life QCD such a scale is given by A x = 47r/ 7r ~ l.lGeV(> Aq CD ), 
above which (Q 2 > A 2 ) the QCD gauge coupling a becomes small enough for the 
OPE to be obviously valid. On the other hand, in the large Nf QCD the QCD gauge 
coupling a is almost constant (walking) over the range A 2 < Q 2 < A 2 and is not very 
small (of order 1), where A x = Airf^ ~ 4.7m(< A) in the large Nf QCD [21] (m 
is the dynamical mass of the fermion). Nevertheless we expect that the OPE is still 
operative as was discussed [IS] in the walking/conformal (scale-invariant) technicolor. 
Note that OPE obviously breaks down for Q 2 < m 2 since the gauge coupling blows 
up due to the fact that the fermion which acquired mass m gets decoupled from the 
beta function, while for Q 2 > A 2 OPE is trivially valid (with vanishing anomalous 
dimension) since the gauge coupling is negligibly small due to the strong asymptotic 
freedom (see Fig. [TJ. 

Then the OPE dictates that the correlator Uyy(Q 2 ) — U AA (Q 2 ) behaves as for 
large Q 2 : 

Tl v _ A (Q 2 ) = Tl vv (Q 2 ) - U AA (Q 2 ) ~ 10b. for Q 2 > A 2 ~ (4vr^) 2 , (6.8) 

where 7 m is the anomalous dimension, and {qq)\ x is qq condensation at the scale of 
A x . 6 Note that the large Nf QCD is a walking theory characterized by 7 m ~ 1 jTSJ 
which was confirmed by explicit calculation [21]. On the other hand, in the real-life 
QCD we have 7 m ~ up to logarithm. 7 Eq. ()6.8|) roughly reflects the slowly damping 
behavior of (Il vv (Q 2 )— Tl AA (Q 2 )) / f 2 of the large Nf QCD compared with the real-life 
QCD in Fig. [H] at least in the high energy region Q 2 > A 2 . In fact, 1/Q 2 damping for 
7 m ~ 1 yields the integral for Am 2 ^ uv ^ in Eq. (|fi.7j) a logarithmic divergence which is 
cutoff by A 2 beyond which the integrand changes into the 1/Q 4 damping due to the 
strongly asymptotic freedom there. Then the enhancement by the large anomalous 
dimension 7 m ~ 1 reads roughly 

Am *(uv) i / a 2 / a 2\ 2g 



SrQQ V) 
P 

J 7T 



ln(A 2 /A 2 ) ~ , * = - 2 ln(47rp), (6.9) 



where we have used Eq. (16. with p ~ 15.1, q ~ 4.9. This implies that Am^ UV j//^ 
diverges as inverse square root scaling, when we approach to the critical point. As 
to the infrared contributions to Am 2 ^ uv ^// 2 , we can see from the areas divided by 
the vertical line of Q 2 = (47if n ) 2 in Fig. Elthat, even evaluated with AN um (^ A), 
Am 2 ^)// 2 dominates Am 2 ^// 2 , 70% to 30%, already at a* = 0.90, in sharp 
contrast to the real-life QCD where UV contributions vs IR contributions are opposite 
ratio: 20% to 80%. We expect that as we approach to the critical point a* — > a CT , 
the cutoff A 2 of the integral region in Fig. [f3 grows rapidly even if the damping slope 

6 Making explicit the dimensions, Eq. JO} reads U V ^ A (Q 2 ) - (A x ) 2 • (A x /Q) 4_27m . 

7 In the real-life QCD the OPE in Eq. I|tj.8f) actually has a logarithmic factor a(Q 2 ) x 
(ln(Q 2 /A| CD )) A , with A = 24/(33-2^/) = 8/9 (N f = 3), which however reads (ln(Q 2 / A q C D))~ 1/9 
and yields only a negligible Q 2 dependence. 
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were unchanged, since the figure is drawn in the unit of f 2 which actually vanishes 
exponentially with respect to A 2 . Then Am^ UV )//^ — > Am 2 // 2 near the critical 
point. Hence we conclude an inverse square root scaling for Am 2 // 2 : 

Am 2 1 

—FT ~ / 6 - 10 



This is compared with our numerical data in Fig. where the solid line is a fitting 
function as 

I(a t ) = m x I 2n +n\ (6.11) 




which was obtained by fitting it to the calculated data. Best fitted values for m and 
n are found to be m ~ 0.035 and n ~ —0.82. Here we note that Fig. was given by 
the contributions from the region < Q 2 < Aj^ um instead of < Q 2 < A 2 . We shall 
later discuss additional contributions from A^ um < Q 2 < A 2 which may be estimated 
analytically based on the walking behavior of Uv-a{Q 2 ) ~ ^/Q 2 - I n spite of lacking 
the contributions from A^, um < Q 2 < A 2 , our data nevertheless show an enhanced 
mass of pseudo NG bosons reflecting an inverse square root scaling of the walking 
theory: 

Am 2 T . , 0.07tt 

~fT ~ J («*) ~ 7 (6-12) 

J - 



near the critical point a* = a*(iV/, N c ) —> a CT 
In passing, Eq. ()6.8|) further suggests that 



rcal-lifeQCD 



Q 4 U V _ A (Q 2 ) _ ~ Q 2 U V _ A {Q 2 ) Ar ~ const, for Q 2 > A 2 (6.13) 



largcAfjQCD 



It is amusing to note that the first relation seems to hold numerically all the way down 
to the infrared region Q 2 < A 2 where the OPE as it stands obviously breaks down, 
see Fig. |7J 8 The remnant of the OPE carrying the information of the anomalous 
dimension seems to persist as the difference of power behaviors according to the 
difference of the anomalous dimension even in the IR region. 



6.4 Calculation of S parameter 

In Fig. 03 we plot S for several values of a* in the range of a* £ [0.89 : 1]. There is an 
ambiguity how many data points we use when we calculate the differential coefficient 
of Tlyy{Q 2 ) — Haa(Q 2 ) at Q 2 = 0. We calculate it by using 3, 4, 5, and 7 data points, 
and resultant values of 5 are indicated in Fig. [H] by o, x, *, and •, respectively. The 
value of S is sensitive to the choice of number of data points we use, however, it 

8 The line of the large Nf QCD in Fig. [7| does not yet level off to a constant value up till this 
momentum region. This is consistent with Ref. where the numerical values of 7,„ calculated for 
the values of a* still away from the critical coupling are slightly larger than 1. 
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Figure 7: Comparison between {r^ ' ny-yi(Q 2 )//^ in the large iVy QCD with 

a, = 0.90 (solid line) and (^-^ . Il v _ A {Q 2 )/f% in real-life QCD with t F = 0.04 
calculated in Ref. (dashed line). 



converges as we increase the number of data points. From this figure, we see that S 
in the large Nf QCD near the phase transition point takes values 



S =0.25 to 0.30, 



(6.14) 



which are slightly smaller (but without much difference) than that of the real-life 
QCD calculated in the same method, S ~ 0.33, and have a tendency decreasing as 
we approach the critical point. 

Although our results are still larger than the estimation from the perturbative 



calculation 



pert 



6-7T 



0.16, at least in the range of a* G [0.89 : 1] investigated 
here, a simple-minded extrapolation seems to suggest that it might get to smaller 
values in the range of a* G [a cr : 0.89]. Although in the present analysis we were not 
able to solve the IBS equation in the range of a* G [a cr : 0.89] because the numerical 
calculation near the critical point is quite difficult, it will motivate future work to 
extend the analysis to the region closer to the critical point and see whether or not 
S gets dramatic reduction near the critical point. 

In passing, here we note an interesting correlation between S and the numerical 
value of AmJ with an artificial cutoff Anudi- In Fig. El AmJ/J^ is proportional to 
the area surrounded by the vertical axis, horizontal axis, the curve of (JlyviQ 2 ) — 
^aa{Q 2 )) I f% and the line of the cutoff Q 2 = Aj5j um which does depend on the slope 
of the current correlators Uv-a(Q 2 ) at Q 2 = 0, namely the S. When S gets smaller, 
A^ um defined by Eq. ()6.3|) necessarily increases and so does Am, / f 2 in our numerical 
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Figure 8: Values of S calculated from Eq. (|3.5(l for several values of a* in the range of 
a* € [0.89 : 1]. Points indicated by o, x, *, and • are values of S calculated by using 
3, 4, 5, and 7 data points of Tlyv — ^aa, respectively. Solid line shows the form of 
function J(a*) in Eq. (|6.16() . Value of S of the real- life QCD calculated in Ref. (23| is 
also shown in the figure. Here, we used the data with the infrared cutoff parameter of 
the running coupling tp = 0.04. 
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Figure 9: Values of S ■ Am 2 // 2 for several values of a*. Solid line shows the value 
obtained by averaging the values of data below a* = 0.95 as shown in Eq. (|6.15|) . Value 
of 5 • Am 2 // 2 of the real-life QCD calculated in Ref. [S3] is also shown in the figure. 
Here, we used the data with the infrared cutoff parameter of the running coupling 
t F = 0.04. 



calculations. 9 In Fig. El we plotted S ■ Ami/ f° r several values of a* in the range 
of a* G [0.89 : 1]. This product is almost constant near the critical point and seems 
to stay as a constant also in the region a* G [a cr : 0.89]: 

Arn^ 

xS~ constant ~ 0.14 (a* < 0.95) (6.15) 

fl 

From this we might infer the behavior of S as an inverse function of Am\j given 
in Eq. (Jfi~TT|l : 

S~ J( a ,) = ^li = 4| 271 - 0.82 ) . (6.16) 

However, this argument is totally based on the assumption that A^ um /A^ scales 
in the same way as A 2 /A^ does in a Miransky scaling. Although there is no solid 
theoretical support for such a continuation of the Miransky scaling in the more vicinity 
of the critical point as we discuss later, A^, um /A 2 imitates the Miransky scaling as far 
as the region we explicitly calculated is concerned. Thus, it is certainly interesting 
to check in future the above square root scaling of S by directly computing S in the 
region closer to the critical point. 



9 This is contrasted to the case where our integration is cut off at A 2 (instead of A^ um ) which has 
no correlation with the slope at Q 2 = 0. 
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In Ref. [2E]> it was argued that there arises a negative contribution proportional 
to the dynamical mass of the fermion due to the walking behavior of the running 
coupling, and that the S parameter in the large Nf QCD could be much reduced 
relative to that in QCD-like theories. However, as far as we analyzed in the present 
paper S does not seem to take a negative value. 



7 Conclusion and Discussions 

In this paper, in the framework of the SD and the IBS equations, we calculated f n , S, 
and AmJ on the same footing in the large Nf QCD, through the difference between 
the vector current correlator Hyv an d the axial-vector current correlator 11^. 

When the chiral phase transition point is approached from the broken phase, 
a* — > a cr , f n and AmJ go to zero both with the essential-singularity scaling (Miransky 
scaling), whereas the ratio AmJ//^ increases and actually scales as an inverse square 
root (Fig. El and Eq. (jgJ2D ): 

Am| 0.07. 

Thus the radiative mass of the pseudo NG boson can be dramatically enhanced com- 
pared with the conventional estimate of the techni pseudo's based on a simple scale-up 
of QCD with a factor scaling of Nf and N c . 

Here we should mention that our numerical estimate of AmJ//^ was done on the 
integral for the region < Q 2 < A^ um instead of < Q 2 < A 2 , since it is rather 
difficult to estimate them by directly solving the IBS equation numerically including 
larger Q 2 regions. We now discuss the additional contributions from A^, um < Q 2 < 
A 2 through an analytical estimate, instead of directly solving the IBS equation, based 
on the observed walking behavior of the current correlators Uv-a{Q 2 ) ~ 1/Q 2 - We 
divide the integral for Am 2 as 

J 7T \ Num / 

4(Num) + ^ m ^(UV>)' 

By definition of Anuhi in Eq. (jfj.Hjl we have 



A ^(Num) + A ^(UV»- ( 7 - 3 ) 



rWQ 2 ) i a 2 

ft 50 Q 2 



V ) " ( A Lm < Q 2 < A 2 ) , (7.4) 



which is substituted into the above Am 2 ( UV> - ) : 

Am Vv>) J_ 3a em Agj um 
ft " 50 4vr ft 




(7.5) 



Values of (A^ um /A 2 ) and (A 2 /A^ um ) are summarized in Table [TJ Considering 
the definition of A x (= Anf^) and the fact that the scaling of f n /A is identified 



22 



a* 


A 2 /A 2 

Num/ x 


A 2 /A2 


a* 


A 2 /A 2 


A 2 /A 2 

I *- / J1 Num 


0.890 


28 


3.8 x 10 5 


0.950 


21 


2.7 x 10 3 


0.895 


29 


2.0 x 10 5 


0.955 


21 


1.9 x 10 3 


0.900 


26 


1.3 x 10 5 


0.960 


21 


1.4 x 10 3 


0.905 


27 


7.6 x 10 4 


0.965 


21 


1.1 x 10 3 


0.910 


26 


4.8 x 10 4 


0.970 


20 


8.8 x 10 2 


0.915 


25 


3.3 x 10 4 


0.975 


20 


6.7 x 10 2 


0.920 


26 


2.0 x 10 4 


0.980 


20 


5.2 x 10 2 


0.925 


25 


1.4 x 10 4 


0.985 


19 


4.1 x 10 2 


0.930 


23 


1.0 x 10 4 


0.990 


19 


3.3 x 10 2 


0.935 


23 


7.0 x 10 3 


0.995 


18 


2.7 x 10 2 


0.940 


23 


4.9 x 10 3 


1.000 


18 


2.2 x 10 2 


0.945 


22 


3.6 x 10 3 









Table 1: Values of (A^ um /A^) and (A 2 /A^ um ) for several values of a*. 



as Miransky scaling (see Eq. ()6.9j) and Fig. EJ), we can understand that A 2 /A 2 also 
shows the (inverse) Miransky scaling. A 2 /A 2 is viewed as a product of A 2 /A^ um and 
A^, um /A 2 , the former becoming more prominent as we get closer to the critical point. 
This implies that cutting off the integral for Am 2 at A.^ um does not reproduce main 
part of the walking theory and hence does not yield a good approximation, in sharp 
contrast to the case of the real-life QCD. Nevertheless A^^/A 2 is also growing as we 

get closer to the critical point, which corresponds to decreasing of S. This actually 
imitates the Miransky scaling in such a way that Am^ Num j behaves as if in an inverse 
square root scaling observed in Fig. El and Eq. ()6.12|) . 

In Fig. [TUl we plot the values of Am^ Num ), Am 2 ( - uv> - ) , and Am 2 in the unit of f 2 
for several values of a*. As was noted above, the additional contributions from A^ um < 
Q 2 < A 2 is much larger than those from Q 2 < A^ um which we calculated numerically, 
since A 2 /A 2 xJum A^ um /A 2 (see Table [TJ in the walking regions. The fitting function 
for total Am 2 reads m ~ 0.042 and n ~ —0.268 in the parameterization of Eq. (|6.11J1 . 
This figure clearly shows that the contribution from Q 2 > A^ um dominates the scaling 
of total Am 2 . When it is applied for a popular technicolor model (one-family model) 
of Farhi-Susskind ^2], with A x = 4nf n ~ An x 123GeV, we are interested in the region 
A < Aetc ~ 10 3 TeV and hence A/A x ~ 6.5 x 10 2 (we are using N TC = 3 in this 
case). This roughly corresponds to a* = 0.975 in Table [T] and Am 2 // 2 = 1.12, about 
9 times larger than the value estimated in the real-life QCD Am 2 // 2 ~ 0.123 |2*3~] 
(experimental value Am 2 = 0.148 ± 0.001) in the same method. This suggests that 
the radiative mass of the pseudo NG bosons in the walking technicolor could be 
enhanced by a factor 3 compared with a simple scale up of the QCD. The typical 
estimate of lOOGeV mass range [12] would be boosted to 300GeV range, which could 
be of phenomenological relevance. 

As to S, our numerical result showed S ~ 0.25 — 0.30 which is somewhat smaller 
than that in the real-life QCD and indicated a decreasing tendency as we approach the 
critical point. Here we note that there is a subtlety to evaluate the slope of YIv-a{Q 2 ) 
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Figure 10: Values of AmL Num j, Ato^ uv> ^, and Ato^ in the unit of for several 
values of a*. 



at Q 2 = numerically. To check this subtlety we compute another quantity S*: 

U V _ A (Q 2 = f 2 ) - Tl v _ A (Q 2 = 0) 



5* 



-47T 



P 

J 7T 



(7.6) 



which is similar to the popular definition of S except for the point that ours is defined 
for the space-like momentum instead of the time-like one. Numerical results for several 
values of a* are summarized in Table [21 

In Fig. ^2 we plotted the values of shown in Tabled The data in the previous 
definition are also plotted in this figure. The result is somewhat larger than the value 
in the previous method and seems to indicate slower decreasing behavior than the 
previous one as we approach the critical point. Thus, up to this subtlety, it seems 
unlikely that S is dramatically reduced in the large Nf QCD at least in the region we 
calculated. At any rate it is highly desirable to extend our analysis to closer region 
to the critical point and check by explicit computation whether or not our conclusion 
in this paper persists there. 

In the present analysis, we only considered SU(N C ) gauge theory with N c = 3 
(with Nf flavors). This is because, by taking N c = 3, it became possible to compare 
the result of the large Nf SU(N C ) gauge theory with usual QCD (iV c = 3, Nf = 3), 
and see the difference of dynamics which is caused by the difference of the number 
of flavors. However, there is no reason to prohibit us from considering the large 
Nf SU(N C ) gauge theory with iV c ^ 3 as an underlying theory of the electroweak 
symmetry breaking. (As an example of recent work on walking technicolor model 
based on SU(2) gauge theory, see So it is interesting to investigate S and 
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a* 


U V . A (Q 2 = 0)/A 2 = f 2 /A 2 


n v - A (Q 2 = f*)/A 2 




0.890 


6.05 x 10~ 1U 


5.90 x 10~ 10 


0.308 


0.895 


1.08 x 10~ 9 


1.05 x 10~ 9 


0.309 


0.900 


1.85 x 10~ 9 


1.80 x 10~ 9 


0.311 


0.905 


3.09 x 10~ 9 


3.01 x 10~ 9 


0.311 


0.910 


5.00 x 10~ 9 


4.88 x 10~ 9 


0.311 


0.915 


7.91 x 10" 9 


7.71 x 10~ 9 


0.312 


0.920 


1.22 x 10~ 8 


1.19 x 10~ 8 


0.313 


0.925 


1.85 x 10~ 8 


1.80 x 10~ 8 


0.314 


0.930 


2.74 x 10~ 8 


2.67 x 10~ 8 


0.314 


0.935 


3.99 x 10~ 8 


3.89 x 10~ 8 


0.315 


0.940 


5.72 x 10~ 8 


5.57 x 10~ 8 


0.316 


0.945 


8.06 x 10~ 8 


7.86 x 10~ 8 


0.316 


0.950 


1.12 x 10~ 7 


1.09 x 10~ 7 


0.317 


0.955 


1.54 x 10~ 7 


1.50 x 10~ 7 


0.318 


0.960 


2.08 x 10" 7 


2.03 x 10~ 7 


0.318 


0.965 


2.78 x 10~ 7 


2.71 x 10~ 7 


0.319 


0.970 


3.68 x 10~ 7 


3.59 x 10~ 7 


0.320 


0.975 


4.82 x 10- y 


4.70 x 10- Y 


0.320 


0.980 


6.25 x 10~ 7 


6.09 x 10~ 7 


0.321 


0.985 


8.02 x 10~ 7 


7.82 x 10~ 7 


0.321 


0.990 


1.02 x lO" 6 


9.96 x 10~ 7 


0.322 


0.995 


1.29 x 10- 6 


1.26 x 10~ 6 


0.322 


1.000 


1.62 x 10- 6 


1.57 x 10~ 6 


0.323 



Table 2: Values of U V -a{Q 2 = 0)/A 2 = f 2 /A 2 , T\ V - A {Q 2 = f 2 )/^ 2 , and & for 
several values of a* . 



Am 2 , in the case of N c ^ 3 in future works. 
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Figure 11: Values of 5 1 * shown in Tabled Values of S are also plotted in the figure. 



iS F (p) 




full fermion propagator 



Figure 12: A graphical expression of the SD equation in the (improved) ladder ap- 
proximation. 

A Appendix : SD equation in the improved ladder 
approximation 

Schwinger-Dyson (SD) equation is a powerful tool to study the dynamical generation 
of the fermion mass directly from QCD (for reviews, see, e.g., Refs. (22 EH])- The 
SD equation for the full fermion propagator iSp 1 = A{p 2 )$ — B(p 2 ) in the improved 
ladder approximation [>35] is given by (see Fig. ^] for a graphical expression) 

(p- q)n{p- q)v 

(p-q) 2 



iS F \p)-i> = C 2 j^ifiPA) 



(p - 



(A.l) 

v is the second casimir invariant, g{p,q) is the running cou- 
pling, and the Landau gauge is adopted for the gauge boson propagator. The SD 
equation provides coupled equations for two functions A and B in the full fermion 
propagator Sf- When we adopt a simple ansatz for the running coupling, g 2 (p,q) = 
<7 2 (max(p|5, q%)) [33], with (p 2 E ,q E ) being the Euclidean momenta, we can carry out 



where C 2 (- 



26 



the angular integration and get A{p 2 ) = 1 in the Landau gauge. Then the SD equa- 
tion becomes a self-consistent equation for the mass function S(p 2 ) = B(p 2 ). The 
resultant asymptotic behavior of the dynamical mass S(p 2 ) is shown to coincide with 
that obtained by the operator product expansion technique. 

However, it was shown in Ref. [12] that the axial Ward-Takahashi identity is 
violated in the improved ladder approximation unless the gluon momentum is used 
as the argument of the running coupling as g 2 {(pE ~ Qe) 2 )- in this choice we cannot 
carry out the angle integration analytically since the running coupling depends on 
the angle factor cos9 = pe ■ Qe/\pe\\qe\- Furthermore, we would need to introduce a 
nonlocal gauge fixing [1011121 to preserve the condition A = 1. 

In Ref. (20], however, it was shown that an angle averaged form g 2 (p E + q E ) gives 
a good approximation. Then, in the present analysis we take the argument of the 
running coupling as 

9 2 (PE,q E ) => fipl + ql)- (A.2) 

After applying this angle average approximation and carrying out the angular inte- 
gration, we can show (see, e.g., Refs. [IE!) that A always satisfies A{p 2 ) = 1 in the 
Landau gauge. Then the SD equation becomes 

Ws) = C 2 I dy V S ^ ^ + V) ( A 3 ) 

lQn 2 J y + T l 2 (y) max(x,y) ' 

where x = p\ and y — q% . 



B Appendix : Numerical method for solving the 
SD equation 

In this appendix we briefly explain how to solve the SD equation numerically. 

We first introduce the infrared (IR) cutoff Xsd and ultraviolet (UV) cutoff A$d 

as 

A 2 e XsD/A < x, y < A 2 e As ° /A . (B.l) 
Then, we discretize the momentum variable x and y into Nsd points as 

x, = A 2 exp [A SD /A + D SD • i ] , (i = 0, 1, 2, • • ■ , (N SD - 1)) , (B.2) 

where 

„ _ (Asp - A so )/A 

~ N SD -! ■ (&3) 
Accordingly, the integration over y is replaced with a summation as 

dy => Dsd^Vj ■ (B.4) 



3 

Then, the SD equation in Eq. ()A.3|) with the running coupling in Eq. (|A.2|) is rewritten 

as 

EW = D » ? ^ + ^b~) ^k)- (R5) 



27 



This discretized version of the SD equation is solved by the recursion relation: 

*W**) = -V, 7) »V, • //,! P—, , ■ (B.6) 

^ 2 j max(xi, yj) yj + £(„)(%•) 

Starting from a suitable initial condition (we choose S( )(xj) = 1), we update the 
mass function by the above recursion relation. Then, we stop the iteration when the 
convergence condition 



2 



2 



is satisfied for sufficiently small e, and regard this £(„) as a solution of Eq. (jB.5|) . 

In the present paper, we use the following parameters to solve the SD equation 
numerically: 

A 5D /A = +24 , A 5 d/A = -24 , N SD = 1000 , e = HT 20 . (B.8) 



C Appendix : Numerical method for solving the 
IBS equation 

In this appendix we transform the IBS equation in Eq. (j5.1|) into the form in which 
we can solve it numerically. 

First, we introduce the conjugate bispinor bases defined by 



l1 J) 0; q, e) = 7 r^ j (p*; q, e) f 7 



(J), 



(C.l) 



Multiplying these conjugate bispinor bases from left, taking the trace of spinor indices 
and summing over the polarizations, we rewrite Eq. (|5.1|) into the following form: 



TV\u,x)x?\u,x) 



1X' J J-oo 



dv j dyy 2 K^\u,x;v,y)x ( j J \v,y) = 4 J) (w,x), 

(C.2) 



where the summation over the index j is understood, and 



J> [Ti J \p;q,e) (e ■ G {J) 



T>p(u,x 



) = ^tr[rWfaq,e)Ttoq)r<% ]q ,e) 



K-p{u,x;v,y) = I dcos6 ^ tr [t\ J \p; q, e)K(p, k)Tf\k; q, e) 
J— l e 

with the real variables v and y introduced as 

k ■ q = —v Q, k ■ p = —uv — xy cos 6. 
Here 9 is the angle between the spatial components of p^ and k^. 



(C.3) 
(C.4) 

(C.5) 
(C.6) 
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Using the property of xY^ i n Eq. (|4.1()jl . we restrict the integration range as v > 0: 



dvK lj (u,x;v,y)Xj J \v,y) = / dv [K^u, x; v, y) + A%(u, x; -v, y)) Xj J \v, y). 

(C.7) 

Then, in the following, we treat all the variables u, x, v and y as positive values. 

To discretize the variables u, x, v and y we introduce new variables U, X, V and 
Y as 



(C.8) 



u = e , x 
v 

v — e , y = e 

and set ultraviolet (UV) and infrared (IR) cutoffs as 

U,V e [Aa, At/], X,F G [Ax, Ax]. (C.9) 

We discretize the variables {7 and V into Nbs,u points evenly, and X and Y into 
Nbs,x points. Then, the original variables are labeled as 

U[i v \ = exp [At/ + I>t//t/] , X[/ x ] = exp [Ax + D X I X ] , 
V[i v ] = exp [At/ + Dul v ] , = exp [Ax + D x h] , 

where I v , I v = 0, 1, 2, • • • (ATr^ - 1) and 7 X , iy = 0, 1, 2, • • • (N B s,x ~ 1)- The mea- 
sures Dt/ and D x are defined as 

A// — At/ ^ Ax — Ax 



D 



u 



D 



x 



Nbs,u — 1 Nbs,x — 1 

As a result, the integration is converted into the summation: 



/ y 2 dy dv 

Jv>0 



D V D V vy 3 

IvJy 



(CIO) 



(c.n; 



In order to avoid integrable singularities in the kernel K(u, x; v, y) at (u, x) = (v, y), 
we adopt the following four-splitting prescription [ST] : 



K i:j (u,x,v,y) 



- [ Kij(u,x,v + ,y + ) + K ij (u,x,v + ,y-) 
+ Kij(u, x, v-,y + ) + Kij(u, x, v_,y_) ], (C.12) 



where 



v± = exp 



V± 



D 



u 



V± = exp 



Y± 



D 



x 



(C.13) 



Now that all the variables have become discrete and the original integral equation 
(|5.1|) has turned into a linear algebraic one, we are able to deal it numerically. 

When we solve the IBS equation numerically in this paper, we use the following 
parameters : 

[At/, At/] = [ 34.5 x a* -45.0 , 0.0], 
[Ax, Ax] = [ 43.5 x a, -52.0, 0.0], 

N BS ,u = N BS ,x = 34. (C.14) 

Note that At/] and [Ax, Ax] are chosen so that the dominant supports of liyviQ 2 ) — 
^-aa{Q 2 ) always lie within the energy region between UV and IR cutoffs. 
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